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ABSTRACT 

We present results from the fitting of infrared (IR) spectral energy distributions of 21 active galactic 
nuclei (AGN) with clumpy torus models. We compiled high spatial resolution 0.3-0.7 arcsec) 

mid-IR TV-band spectroscopy, Q-band imaging and nuclear near- and mid-IR photometry from the 
literature. Combining these nuclear near- and mid-IR observations, far-IR photometry and clumpy 
torus models, enables us to put constraints on the torus properties and geometry. We divide the 
sample into three types according to the broad line region (BLR) properties; type-Is, type-2s with 
scattered or hidden broad line region (HBLR) previously observed, and type-2s without any published 
HBLR signature (NHBLR). Comparing the torus model parameters gives us the first quantitative 
torus geometrical view for each subgroup. We find that NHBLR AGN have smaller torus opening 
angles and larger covering factors than those of HBLR AGN. This suggests that the chance to observe 
scattered (polarized) flux from the BLR in NHBLR could be reduced by the dual effects of (a) less 
scattering medium due to the reduced scattering volume given the small torus opening angle and (b) 
the increased torus obscuration between the observer and the scattering region. These effects give a 
reasonable explanation for the lack of observed HBLR in some type-2 AGN. 

Subject headings: galaxies: active — galaxies: nuclei — infrared: galaxies 


1. INTRODUCTION 

While active galactic nuclei (AGN) present a variety of 
observational characteristics, the unified model for AGN 
proposes the ubiquitous presence of an obscuring torus 
around their central e ngines, a nd that all AGN are fun¬ 
damentally the same (Antonucci 1993). This optically 
and geometrically thick torus produces the effect of a 
line of sight viewing angle dependency. Type-1 AGN 
are observed with the direct view of fast moving mate¬ 
rial close to the supermassive black hole (SMBH), re¬ 
sulting in broad emission lines in their spectra, while 
type-2 AGN are observed from an edge-on view and the 
torus blocks the broad emission line region (BLR) com¬ 
ponent from our line of sight. The most compelling ev¬ 
idence of the unified model was the detection of polar¬ 
ized broa d em i ssion lines (PBLs) in type-2 AGN (e.g., 


Antonucci & Miller 


1985). Further evidence support¬ 


ing the unified model comes from infrared (IR) obser- 
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vations of several type-2 AGN showing the existence of 
obscured/hidden broad line regions (HBLRs) detectable 
only with dust penetrating infrared observations (e._g., 


Blanco et al. 

1990; 

Nagar et al. |[20021 

2003; Ramos Almeic 

ia et al. 2008] 

)• 


teunanen et al. 


port the unified model, there is the question of why some, 


signs of PBLs. 

Tran 

(2001,2003 

) and 

Moran et al. 

(2001) 

found that only 30 

-507o of type-2 

AGN show 

PBLs. 


Some studies have advocated that the non-detection of 


a PBL i s due to genuine lack of a BLR (e.g., Tran et 
aTpOll). Others have suggested that the non-detection 


is due to obscuration effects, rendering the detection of 
PBLs as difficult or impossible, even with deep near 


IR (N IR) spectro-polarimetric observations (Alexander 
2001| ). Using a statisticall y complete IRAS 61 ) /im se¬ 


lected type-2 AGN catalog, Heisler et al. ([1997) investi¬ 
gated the relationship between the detectability of PBLs 
and IR color as an indicator of the torus inclination angle. 
They showed that only AGN with a low torus inclination 
angle have high detection rate of PBLs compared to those 
with high inclinations. This result strongly suggests that 
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PBLs could be obscured when there is an edge-on view 
through the torus and/or nuclear obscuration in the host 
galaxies. In addition to the optical spectro-polarimetry, 
X-ray observations suggest that there is a weak evidence 
s howin g different absorption in two types of type-2 AGN. 
Gu et al. (2001) found that the AGN with PBL have 
slightly lower c olumn density (AM) th an those without 
PBL. Similarly, Lumsden et al. (2004) showed that the 
detection rate of PBL decreases as a function of AM, sug¬ 
gesting that the absorption effect by dusty torus could 
play a role on the detectability of PBL in AGN. 

To understand the role of the obscuration by the torus 
in type-2 AGN and the detectability of PBLs, knowing 
the torus geometry and properties is crucial. In recent 
years much progress has been made toward understand¬ 
ing the torus geometrical structure. Thanks to the im¬ 
provement of computing power, more physically realis¬ 
tic torus models assuming “clumpy” distributions (called 
clumpy torus models) have been coded by several authors 


^__s) have been co ded by s everal authors 

(iNenkova et al.||2002l |2008a|b[ iHonig et al. 1120061 ISchart- 

mann et M.||Mj8l l llonlg & Kishimoto 20 T U lfet a 


al.||2008| |Honig &; kishimoto||2(JlU(~[5l 

2lrr These models readily reproduce hi, 


et al ||2012|). These models readily reproduce high spa¬ 
tial resolution nuclear NIR to mid-IR (MIR) spectral en- 


compact torus of < 10 pc radius (e.g., 

Ramos Almeida 

et al. 

2009, 2011a; Nikutta et al. 2009 

Alonso-Herrero 

et al. 

20111 Honig et al. 

201H Lira et al. 2013). On 

the other hand, traditionaJ 

smooth torus models (Pier & 

Krolik 19921 1993 Efstathiou & Rowan-Robinson|1995| 


of nearby AGN (e.g., Alonso-Herrero et all 2003 Gandhi 


Asmus et al.||2(Jll| Ichikawa et al. |2012ap 
e torus morphology remains far from beim 


et al. 2009 

Still, the true torus morphology' remains far from being 
conclusively determined until future observations resolve 
the torus. Both smooth and clumpy torus models have 
degeneracies, and from SED-fit ting alo ne is not possible 
to choose one or the other (see Feltre et al. 2012, for a 
comparison between smooth anaclumpy torus models). 

In this paper, under the assumption that the torus fol¬ 
lows a clumpy distribution of dust, which we consider 
more realistic in principle , we discuss how the precise 
modeled torus morphology plays a key role in the prob¬ 
ability of the detection of PBL by fitting clumpy torus 
models to our series of IR SEDs. We use 21 high spa¬ 
tial resolution MIR spectra in combination with NIR to 
far-IR (FIR) photometry, constituting one of the largest 
compilations of nuclear IR SEDs of AGN in the lo¬ 
cal universe. These SEDs afforded by 8 m class tele¬ 
scopes minimize contamination of the MIR torus spec¬ 
tra from surrounded diffuse MIR emission from warm 
dust and/or starformation region emission, and hence 
we are able to construct the highest fidelity torus SED. 
Therefore we can probe the detectability of the PBL with 
the least amount of host galaxy contaminations than 
ever before. The luminosity dependence of the torus 
morphology will be discussed in a forthcoming paper 
(Ichikawa et al. in prep.). 

2. OBSERVATIONS 
2.1. The Sample 

Our principal motivation in this study is to investigate 
whether the torus model morphology plays a major role 
in the chance of PBL detections. To achieve this goal, we 


compiled the nearby AGN sources from the MIR sam ples 
of[Go nzalez-Martm et al.|(|2Q13|) (21 sources) and|Alonso- 
Herrero et al.| Q201 ip (13 sources) as both samples already 
compiled the currently available data set of ground-based 
N band spectroscopy. Of the 34 sources, five are in com¬ 
mon in both samples, therefore the total number is 29 
sources. We further set the criteria for survey inclusion 
that the objects must have at least one high spatial reso¬ 
lution NIR (1-5 /im) measurement, as the NIR ba nds sig¬ 
nificantly help to con strain the torus parameters (Ramos 
Almeida et al. 2014). 22 out of 29 sources fulfilled this 
criterion. We also removed NGC 1808 from this study as 
controversy remains to whether it hosts AGN or ultra- 
luminous X-ray sources in the galactic center due to the 
low X-ray luminosity logL 2 -in kev = 40.4 erg s -1 (|Scar- 


rott et al. 1993 Jimenez-Bailon et al. 2005). We sum¬ 
marize the properties of the 21 sources m Table 1. 

Our sample spans AGN bolometric luminosities taken 

from the literature ( L j^) in the range log =42.7- 


45.1 erg s -1 (see Table 1), with a mean value of 
44.0 erg s _1 . This value is fairly cons iste nt with that 
of magnitude-limited S eyfert catalogs (Maiolino & Rieke 


1995 


Ho et al. 


1997). This suggests that our sample 


could be representative of AGN and their tori in the lo¬ 
cal universe, although the sample is not complete. 

2.2. New Observations 

We obtained N (Si2 filter; the central wavelength with 
A c = 8.73 /am and 50% cutoff range of AA = 0.39 /im) 
and Q (Qa filter; A c = 18.06 /im and AA = 0.76 /im) 
band imaging data of NGC 5135 and NGC 5643, ob¬ 
served by T-ReCS (Program ID GS-2012A-Q-43, PI: 
Nancy Levenson). The standard MIR chop-nod tech¬ 
nique was performed for the observations. The data re- 
duction wa s made by using RedCan (Gonzalez-Martfn 
etaLp0l3 ). 


2.3. Published Data from the Literature 

We collected the estimated values of the nuclear NIR to 
MIR emission when available. We compiled nuclear NIR 
data from both ground- and space-based telescopes such 
as VLT/NACO and iKT/NICMOS. The only exception 
is the galaxy NGC 5728, whose only NIR flux is from 
2MASS ( |Peng et al.||2006~| ), which we use as an upper 
limit. Information about the NIR fluxes used here is 
compiled in Table 2 (columns 2-6). 

In order to further reduce the parameters space at 
longer wavelengths, we used the Spitzer /'IRS 30 /im con¬ 
tinuum fluxes reported by Deo et al. (2009) as upper 


limits in our fits. We do that because the star forma¬ 
tion component can be important at 20-30 /im, and 
beyond 30 /im complet ely overwhelm s the AGN torus 
emission in most cases (Netzer et al. 2007) at the spa¬ 
tial resolutions afforded by Spitzer. We only consider 
the Spitzer 30 /im photometry as data points, and also 
included the Spitzer /IRS spectroscopy when the AGN 
spectral turnover at 20-30 /im is clearly seen in the IRS 
spectra. This feature suggests that the torus emission is 
dominant even in the large aperture data from Spitzer 
(Alonso-Herrero et al. 2012b Honig et al. 2014). Only 
two sources fulfilled the criterion (1C 4329A and MCG - 
5-23-16). For those galaxies, we also collected available 
Herschel /PACS data. All the FIR flux information is 
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TABLE 1 

Properties of the Sample 


Name 

z 

d 

Slit/Size 

Type 

Group 

TVh 

log 4^ 

b/a 

Ay 

i 

Ref 

(i) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

(9) 

(10) 

(ii) 

(12) 

NGC 1365 

0.0055 

18 

0.35/31 

Syl.8 

Type-1 

23.6 

42.9 

0.5 

<5 


(A1,B1,B1,- • •) 

NGC 4151 

0.0033 

13 

0.36/23 

Syl.5 

Type-1 

22.8 

43.7 

0.71 



(A9,A9, •••,•••) 

IC 4329A 

0.016 

65 

0.75/240 

Syl.2 

Type-1 

21.8 

43.6 

0.28 



(A10,A9,•••,•••) 

NGC 7469 

0.016 

66 

0.75/240 

Syl 

Type-1 

20.7 

45.1 

0.72 



(A9,A9,- •■,•••) 

NGC 1068 

0.0038 

15 

0.36/26 

Sy2 

HBLR 

> 25 

45.0 

0.85 


60-90 

(A2,A9,- ■ • ,A9) 

NGC 2110 

0.0078 

31 

0.36/54 

Sy2 

HBLR 

22.5 

43.9 

0.74 

5 

40 

(A9,A9,A9,A9) 

MCG -5-23-16 

0.0085 

34 

0.75/120 

Sy2 

HBLR 

22.2 

44.4 

0.46 

>6 

53 

(A9,A9,A9,A9) 

NGC 3081 

0.008 

32 

0.65/100 

Sy2 

HBLR 

23.9 

43.8 

0.8 



(A3,B2, •••,•••) 

NGC 3227 

0.0039 

17 

0.75/62 

Sy2 

HBLR 

22.2 

43.4 

0.68 



(A11.A9, •••,•••) 
(A8,A9,A9,A9) 

Circinus 

0.0014 

4 

0.60/12 

Sy2 

HBLR 

24.6 

43.6 

0.44 

9 

60-90 

NGC 5506 

0.0062 

25 

0.36/44 

Sy2 

HBLR 

22.4 

44.2 

0.30 

>11 

40 

(A9,A9,A9,A9) 

IC 5063 

0.011 

46 

0.67/150 

Sy2 

HBLR 

23.3 

44.5 

0.68 

7 


(A2,A9,A9,- • •) 

NGC 7582 

0.0053 

21 

0.75/76 

Sy2 

HBLR 

22.7 

43.3 

0.42 

8,13 


(A9,A9,A9,- ••) 

NGC 7674 

0.029 

118 

0.75/430 

Sy2 

HBLR 

>25 

45.0 

0.91 

-3-5 


(A9,A9,A9,- ••) 

NGC 1386 

0.0029 

11 

0.31/17 

Sy2 

NHBLR 

>25.0 

42.9 

0.4 


65,85 

(A2,B2,- • ,C1) 

NGC 3281 

0.011 

43 

0.35/73 

Sy2 

NHBLR 

24.3 

44.6 

0.4 



(A4,B1,•••,•••) 

Cen A 

0.0018 

3 

0.65/11 

Sy2 

NHBLR 

23.7 

44.0 

0.4 

14.0 


(A5,B2,A9,- ■ •) 

NGC 5135 

0.014 

59 

0.70/200 

Sy2 

NHBLR 

>25.0 

44.4 

0.7 



(A2,B2, 

NGC 5643 

0.004 

16 

0.35/29 

Sy2 

NHBLR 

23.8 

42.7 

0.9 



(A6,B5,•••,•••) 

NGC 5728 

0.0094 

40 

0.35/69 

Sy2 

NHBLR 

23.6 

44.5 

0.6 



(A7,B6,- ••,•••) 

NGC 7172 

0.0087 

35 

0.36/61 

Sy2 

NHBLR 

22.9 

43.8 

0.46 



(A2,A9,•••,•••) 


Note. — Sample properties. The sample is divided into three subgroups with type-l/HBLR/NHBLR respectively from top to bottom. 
(T) object n am e; (2) redshift; (3 ) luminosity distance (Mpc) gathered from literature for the case of nearby sources. Within the sample 
ofjGonzalez-Martfn et al. (120 13b for NGC 1365. NGC 138 6. NGC 1808, NGC3081, NGC 3281, and Cen A, the values of distance t o the 
galaxies Have been takentrom Ramos Almeida et al.|(|2009). For N GC 5643, the distance has be en taken from Guainazzi et al.|(|2004). For 
the sample of Alonso-Herrero et al. (20 Tl]>, we gathered them from | Alonso-Herrero et ah I pm} . For the other sources, we calculated the 
distances by using cosmological parameter Ho = 75km s _1 Mpc _1 ; (4) slit width (arcsec) / physical size (pc); (5) Seyfert class of AGN. (6) 
Sub group of AGN. Type-1 represents type-1 AGN (Sy 1 to Sy 1.9) based on optical spectroscopy. HBLR represents type-2 (Sy2) AGN with 
hidden broad line region signs, and NHBLR represents type-2 AGN without any published hidden broad line regions signs. (7) hydrogen 
colum n density; (8) logarithm of bolometric luminosit y (erg/s) which is taken from |Gonzalez-Martfn et al.| (|2013|); | Alonso-Herrero et al.| 
([2011). We use a typical bolometric correc tion of 20 (Elvis et al. 1994). (9) The axial ratio; the ratio ot the minor to major axis ot the 
host galaxies. All i nformation is taken fro m Gonzalez-Martin et al. (2U13); (10) Foreground extinction in the unit of mag; (11) inclination 
angle of the torus. Levenson et al. (|2006|) derives the viewing angle ot accretion disk of NGC 1386 and we here assume that the accretion 
disk and the torus are located in the sam e plane ; (12) References of colum n (6b (7b (TO), an d (11) . “• • • ” denotes no r eference. 
References 


et al. 1(11999) 


& Veron 


Itoli et a 



; (B5 


Imanishi (2002); (Bl) Tueller et al. (2008); (B2) Marinucci 
Uuamazzi et al(|( |200T ; (h>6) [Goulding et al. ( 2012]) ; (ui)[l f 


77 al. pOOgfr 


tabulated in column 9 to 11 in Table 2. 

The er ro rs were estimat ed using the prescription given 
by Alonso-Herrero et al. (2012a). For the NACO AO 
observations, we used 20% in J band and 15% in the 
HKLM band. For the other ground-based observation 
data, we applied 30% for J band, 25% for H and K band, 
and 20% for the L band. These errors include the pho¬ 
tometric error, the background subtraction uncertainty, 
and the uncertainty from estimating the unresolved flux. 
M band fluxes were always used as upper limits due to 
the difficulties of estimating the unresolved component. 
For the NICMOS observations, we used 20% for the J 
band, 20% for the H and K band. For the N and Q 
band,we use 15% and 25% errors, respectively. 

2.4. Subsample 

To examine the torus model properties of different 
AGN populations, we divide the sample into subgroups 
based on whether or not the source has HBLR signs in 
previously published observations. We first divide the 
sample into type-1 and type-2 AGN. Althou gh Seyfert 
1.8 /1.9 are very ambiguous objects (e.g., see Elitzur et 


al. 2014 for the details), here we define type-1 as AGIN 


tical spectra. Therefore, we consider Seyfert 1 to 1.9 
as type-1 AGN and Seyfert 2 as type-2 AGN. Next, we 
divide the type-2 AGN into those with any published 
polarized BLR detections in the optical and/or in the 
NIR (HBLR) and t hose w ithout (non-HBLR; hereafter 
NHBLR). We use Marin| (|2014|), who compiled almost 
all the previously published polarization information of 
nearby AGN. These spectro-polarimetric data are taken 
from seve ral large su rve ys in cluding the infrared-selected 
sampl e of|Heisler et al . ([1997 D , the FIR flux limited sam¬ 


ple of A lor an et al. 
optica. 


pie of|Lumsd en et al.FP _ 

(|2000[ |2002p, an d the heterogeneou s 
3lected 


2001), the distance-limited sam- 


and 1V11R selected sample of 


Tran 


Leteroge 

( 2001 ,: 


2003) 


which have at least one broad emission lines in their op- 


Mainly the spectro-polarimetric observations were con¬ 
ducted with small or medium size telescopes (up to 4 m- 
class), while only NGC 3081 has been confirmed to have 
HBLR features with the Keck 10 m telescope ([Moran et 
ahpOOO). Therefore, we should note that some HBLR 
AGN could contaminate the subgroup of NHBLR in 
the cases where the BLR is below the signal-to-noise 
afforded by the 4 m class telescope observations (see 
Ramos Almeida et ah, in prep.). Some sources have cur¬ 
rently no published spectro-polarimetric data, but have 
clear broad emission lines in NIR wavelengths. These 
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TABLE 2 

List of Photometry 


Name 

(i) 

j 

(2) 

H 

(3) 

K 

(4) 

L 

(5) 

M 

(6) 

N 

(7) 

Q 

(8) 

30 fim 
(9) 

70 n m 
(10) 

160 fim 

(ii) 

Ref 

(12) 

Type-1 

NGC 1365 


8.3 ±0.83 




203 ± 30.0 

818 ± 204 

< 12.3 



(A1,B1) 

NGC 4151 

69 ± 14 

104 ± 10.4 

178 ± 17.8 

< 325 


1320 ± 198 

3200 ± 800 

< 3.64 



(A7,A7) 

IC 4329A 


50.0 ±8.0 

102 ± 10 

< 210 


1014 ± 150 


1.52 ±0.015 < 1.79 

< 0.97 

(A7,A7) 

NGC 7469 

16 ±3.2 

40 ±4.0 

68 ±6.8 

< 84 


506 ± 76 

1350 ± 340 




(A7,A7) 

HBLR 

NGC 1068 

9.8 ±2.0 

98.0 ± 15.0 

445 ± 100 

920 ± 140 

2270 ± 340 

10000 ± 1500 

21800 ± 5400 




(A7,A7) 

NGC 2110 




< 33.0 

< 198 

294 ± 44 

561 ± 140 

< 0.8 



(A7,A7) 

MCG -5-23-16 

1.1 ±0.33 

3.7 ±0.93 

10.7 ±2.7 

79.5 ±16.0 

< 139.4 

633 ± 95 

1450 ± 360 


< 1.45 

< 0.45 

(A7,A7) 

NGC 3081 


0.22 ± 0.04 




83 ± 12.5 

231 ± 57.8 

< 1.10 



(A2,B1) 

NGC 3227 


11 ± 1.1 

23 ±2.3 

< 47.0 


320 ± 48 

1100 ±275 

< 1.76 



(A7,A7) 

Circinus 

< 1.60 

4.77 ±0.72 

19 ±2.9 

380 ± 57 

1900 ± 285 

5600 ± 840 

12800 ± 3200 




(A7,A7) 

NGC 5506 

13 ±3.0 

53 ±8.0 

80 ± 12 

290 ± 44 

< 530 

900 ± 135 

2200 ± 550 

< 4.05 



(A7,A7) 

IC 5063 


0.3 ±0.1 

4.8 ± 1.0 



925 ± 139 


< 3.89 



(A7,A7) 

NGC 7582 


11.0 ±1.6 

18.0 ±2.7 

96.0 ±14.4 

141 ±21 

384 ± 57 

527 ± 132 




(A7,A7) 

NGC 7674 

1.25 ±0.25 

5.0 ±0.5 

12.3 ±3.1 

53.0 ±11.0 

< 108 

518 ± 78 


< 1.83 



(A7,A7) 

NHBLR 

NGC 1386 


0.2 ± 0.04 




147 ± 22.1 

457 ± 114 

< 1.58 



(A2,B1) 

NGC 3281 


1.3 ±0.33 

7.7 ± 1.93 

103 ±20.6 

< 207 

355 ± 53.3 

1110 ±278 




(A4,B1) 

Cen A 

1.3 ±0.26 

4.5 ±0.68 

34 ±5.1 

200 ± 30 


710 ± 107 

2630 ± 658 




(A5,B1) 

NGC 5135 

< 0.72 

0.66 ± 0.07 




56.36 ± 8.454 

218.92 ±54.73 

< 3.03 



(A2;A6, B2) 

NGC 5643 


< 1.7 




101.31 ± 15.17 

883.2 ± 220.8 




(A2,B2) 

NGC 5728 



< 7.1 



25 ± 3.75 

184 ± 46 




(A6,B1) 

NGC 7172 


< 0.4 

3.4 ±0.86 

30 ±6 

< 61.4 

165 ± 30 


< 0.98 



(A7,A7) 


Note. — NIR to FIR fluxes used as inputs for BayesClumpy. 
Column 9 corresponds to Spitzer /IRS 30 fi m continuum fl uxes f rom 
10 and 11 list the Herschel /PACS photometry from 


information. The r eferences from NIR to MIR band f 
References . (Al) Carollo et al 


|et al.| |2007| ); (A6) |Heng et al. 


Melendez et a. 
uxes are tabula' 


jym mJy : 

D eo et a ll (12009 b which are used as upper limits in the fits. Columns 
pmfj r and are also used as upper limits. “• • • ” represents no flux 
;ea at column 12. 


(2002) (A2) Quillen et al. (200l|; (A3) Callia no Sz Alloin (20 081); (A4)|Simp son (1998); (A5) iMeisenheimerl 
000| ; (A7) | Alon so-Herrero et al. (2011); (El) Ramos Almeida et al. (2009); (B2) This work. 


sources are M CG-5-23 -16, NGC 2110, and NGC 7582 
(Alonso-Herrero et al.||2011). All the references used for 
dividing the sample Into each subgroup are indicated in 
column 12 in Table 1. 

Finally, the sources in this study are categorized into 
three groups (type-1, HBLR, and NHBLR; see column 6 
in Table 1). The sample contains 4 type-1, 10 HBLR, 
and 7 NHBLR AGN. 


3. APPLICATION OF TORUS MODEL 
3.1. Clumpy Torus Model 


We fit the clumpy torus models of jNenkova et al. 
(2008a), known as CLUMPY , to the data using ~a 
Baye sian approach (B ayesClumpy; Asensio Ramos & 
Ramos Almeida [2009). Here we describe the six free 
CL (JMFY mo del parameters used for the SED fitting 
and the model set-up, which are listed in Table 3. The 
torus clumps are distributed in a smooth, rather than 
sharp, toroidal-shaped boundary of angular width a. 
The inner radius (ri n ) of the torus is set by the location 
of the dust at the sublimation temperature (T su b ^ 1500 
K). This is computed using the AGN bolometric lumi¬ 
nosity L bo i(AGN) 


( ^boi(AGN) 

VlO 45 erg s -1 


0.5 

pc. 


(i) 


The torus has a radial extent (T) defined by Y = 
Amt Am, where r out is the outer radius of the torus. The 
average number of clouds along the line of sight (Alos) 


at a viewing angle i is set as 


Alos = A 0 exp 


" (90-i) 2 ' 


( 2 ) 


where Ao is the average number of clouds along the radial 
equatorial ray. Alos allows us to derive the escape prob¬ 
ability of photons from the AGN central engines (P esc ). 
In the CLUMPY dust distribution, the classification of 
type-1 or type-2 AGN depends on whether or not there 
is a clump along the line of sight, which is a function 
of the viewing angle of the torus, the number of clumps 
and the torus width. This is different from smooth torus 
models, for which the classification of an AGN as type-1 
or type-2 is solely determined by the viewing angle. The 
escape probability of photons passing through the torus 
at a given viewing angle (i) can be calculated as 

Pe sc~e~ NLOS . (3) 


In the CLUMPY model, the radiative transfer equations 
are solved for each clump and thus the calculations de¬ 
pend on the clump distribution within the torus, the 
optical depth of each clump, and also its dust compo¬ 
sition. Here we assume each clump has the same optical 
depth (tv), which is defined at the optical V band. The 
CLUMPY model applies a standard cold oxygen-rich in¬ 


terstellar medium dust, which is called OHMc dust (Os- 


senkopf et al. 1992). The torus clumps are distributed 


as a power law with index q as a function of radius, 
A(r) oc r~ q . 

In addition to these six physical parameters, we add 
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TABLE 3 

Free Parameters of the BayesClumpy 


Parameters 

Parameter range 

Torus radial thickness ( Y ) 

[5, 30] 

Torus angular width ( cr ) 

[15°, 70°] 

Number of clouds along an equatorial ray (Wo) 

[1.15] 

Index of the radial density profile ( q ) 

[0,3] 

Viewing angle ( i ) 

[0°, 90°] 

Optical depth of each cloud (tv) 

[5,150] 


Note. — Torus radial thickness Y is defined as Y = r out /ri n , where 
v o U t is the outer radius and ri n is the inner radius. The cloud distri¬ 
bution between r ou t and ri n is parameterized as r~ q . 


two additional parameters to be fitted or fixed. The 
first parameter is the foreground extinction (Ay), unre¬ 
lated to the torus. Some authors demonstrated that some 
AGN have an extremely deep 9.7 fi m silicate absorption 


obscuration (Alonso-Herrero et al. 

2003| Pohetta et al. 

2008; Alonso-Herrero et al. 

2011; Goulding et al. 2012; 

Gonzalez-Martin et al. 2013 

). They 

■ suggested that dust 


the observed SED and silicate feature absorption. 10 
out of 21 sources are inclined galaxies with low minor- 
to-major axis ratios ( b/a < 0.5; see Table 1). There¬ 
fore, some portion of the observed SED is accounted 


from b y co ol foreground dust extinction. Alonso-Herrero 


et al. (2011) discussed this issue and concluded that for 
Ay > 5, the effects of foreground extinction cannot be 
ignored for reproducing the silicate 9.7 fi m feature. We 
gathered available values of foreground Ay from the lit¬ 
erature and compiled them in column 10 in Table 1. The 
other additional parameter accounts for the multiplica¬ 
tive factor that has to be applied to match the fluxes of 
a given model to an observed SED. Deriving this factor 


enables us to cal 


minosity (Nenkova et al. 


ulate the mod 


Alonso-Herrero et al. 


AGN bolometric lu- 
2008b). As shown by 


(2011), Tboi reproduces well the 
values of L ^ for Seyfert galaxies, and therefore, in the 

following we will refer to L^ d ^ as the bolometric lumi¬ 
nosities of the sample studied here. 

By combining the derived output parameters of the 
CLUMPY model, we can derive other important torus 
morphological parameters as the torus outer radius r ou t, 
defined as: 

Anit = r in Y pc. (4) 

We can also calculate the torus scale height H as: 

H = r ou t sin a pc. (5) 

Finally, we define the “geometrical” torus covering fac¬ 
tor, which is unaffected by the viewing angle, and it is 
defined by integrating the AGN escape probability over 
all angles ( Nenkova et al.|2008a ). This can be written as 


Ct = 1 — 


/•tt/2 

Jo 


PescW) cos(/3)d/3, 


(6) 


where f3 = n/2 — i. Considering that our motivation is 
to characterize the intrinsic torus morphology, the “geo¬ 
metrical” torus covering factor is more relevant here than 
the apparent covering factor. 


3.2. BayesClumpy and Modeling Details 

The CLUMPY database currently contains more than 
5 x 10 6 models. Therefore, when fitting the models to 
the observations, inherent degeneracies have to be taken 
into account. We then use the B ayesClum py fitting 
tool (Asensio Ramos & Ramos Almeida|2009), as it does 
a fast synthesis of the CLUMPY SEDs' in the last ver¬ 
sion of BayesClumpy the inference over the model pa¬ 
rameters can be done either using neural network inter¬ 
polation or multilinear interp olatio n in the full database. 


After running several tests, Ramos Almeida et al. (2014) 


concluded that the latter interpolation produces more ro 
bust results. Therefore, here we use linear interpolation, 
which results in slight differences in the fitted parameters 
(within lcr for the majorit y of t he fi ts) for t he 13 galax¬ 


ies that were modeled by Alonso-Herrero et al. (2011) 
using the neural network interpolation and subsequently 
re-fitted in this paper. 

BayesClumpy can be used to fit photometry and/or 
spectra in a Bayesian scheme, carrying out inference over 
the model parameters for observed SEDs. This way we 
can specify a-priori information about the model param¬ 
eters. Here we consider uniform prior distribution in the 
range of each parameter, as summarized in Table 3. The 
prior distribution of inclination angle (i) is fixed from 
previous observations if avai lable in th e literature, follow¬ 
ing the same approach as in Alonso-Her rero et al. (2011). 
From the ob j ects i n our sample taken from |Gonz alez- 


(2013), NGC 138 6 has two possi ble inclina- 
and 85° (iLevenson et al.|2006). Thus 

an O nno 


Martin et al. 
tion angles 65 
use a uniform prior in the r ange 60 u -90' 


we 

lor this source. 


For the galaxies taken from Alonso-Herrero et al. (2011) 
we use the same inclination angle constraints they em 
ployed, which are compiled in column 11 in Table 1. 

We finally include the AGN direct emission (i.e. a bro 


ken po wer-law ) which is defined in Eq. (13) of Nenkova 
( 2008a| ) to the SED for type-1 AGN, in order to 


et al. 


reproduce the flatte r slope of the NIR band (see Ramos 
Almeida et al.|[2014 for further details). 


4. RESULTS AND DISCUSSIONS 
4.1. Infrared SEDs with BayesClumpy Fitting 

The results of the fitting process to the IR SEDs are 
the posterior distributions for the six parameters that 
describe the model (defined in Table 3), the foreground 
extinction and the multiplicative factor needed to match 
the SED fluxes. However, we can also translate the re¬ 
sults into corresponding spectra, as shown in Figure 1. 

Figure 1 shows the observed SEDs and nuclear MIR 
spectra (black filled dots) of the galaxies with the best 
fit results overlaid (blue solid lines), based on the infer¬ 
ence done with BayesClumpy. The fitted models corre¬ 
spond to those described by the median of the posterior 
distribution of each parameter. All the derived torus pa¬ 
rameters obtained from BayesClumpy are presented in 
Table 4. 

Some SEDs show lower Q-band fluxes than those pre¬ 
dicted by the fitted model. This effect is prominent when 
silicate 9.7 /am feature is observed in deep absorption. 
Although this may suggest that the model spectra still 
have a difficulty to reproduce the 18 /am silicate feature, 
the difference is only within a factor of 3 in the w orst 
case (NGC 3281). See Ramos Almeida et al. (2009) for 
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wavelength [/im] wavelength [/xm] wavelength [/xm] 
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.HBLR. 
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wavelength [/xm] 



wavelength [/xm] 



Fig. 1.— Clumpy torus model fits. The filled dots are the photometric data and the black line is the MIR spectrum. The upper limit 
points are shown as arrows. The solid blue lines are the models computed with the median value of the probability distribution of each 
parameter. The blue shaded areas indicate the range of models compatible with a 68% confidence interval. For the details on the calculation 
of the median values, see Section 3. 
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TABLE 4 

Fitted torus model parameters from SED + spectroscopy data 


Galaxy 

iTtorus 

Y 

N 0 

q 

TV 

i 

C'j' 

1 °g z 'b” od) 

Tin 

T out 

H 


[deg] 





[deg] 


[erg/s] 

[pc] 

[Pc] 

[Pc] 

Type-1 

NGC 1365 

19t£ 

24± 3 

9+ 2 

y -3 

Q O + 0.3 
U,,3 -0.2 

79ls5 

iflli! 

0 16+ 0 - 08 
U -- LO _0.04 

43 -2io.i 

D n^p+°- 005 
U.UOo— 0.005 

1 o + 0-2 
1 -°-0.3 

0-411:1 

NGC 4151 

16+J 

1911 

13±f 

1 g + 0 - 3 
- L-D —0.3 

00 00 

+ 1 
C5 

00 

nil 

0 1 p+ 003 
U - lo -0.02 

43 -9lS'.i 

n 1 nsT 0 • 007 
u.iu»_ 0 008 

9 1 + 0.5 
Z,± -0.5 

o.6ii:i 

IC 4329A 

40+f 

sli 

12±\ 

Q C + 0.1 

148+2 

4+ 4 

^-3 

0 5«+°- 04 
U.OO_o 04 

44 4+ 1 - 8 

0 192+ 0 ' 002 
u.iuz — o 001 

1 - 811:1 

1 i+o-o 
4 - 4 -o.o 

NGC 7469 

2%i 2 2 

to 

to 

1 + 

13±{ 

-i q + 0.3 
1 -' 3 -0.3 

124ll| 

5911 

0 20+°' 05 
U.ZU_o 04 

44 - 6 -°o\ 8 i 

0 239+ 0 ’ 020 

u.zaU-q 021 

5.311.1 

1 9 + 0 ' 4 

4 - y -0.4 

HBLR 

NGC 1068 

56tL 

6±l 

5ll 

0 6+ 1 ' 2 
U.O_o 4 

38ty 

671“ 

0.78+o‘.i9 

44 - 4 ^o.i 

n iqs + 0 - 015 
U. 19 o _ 0 0Q6 

1 0 + 0.6 

4 ■' 4 — 0.1 

1-Uo.l 

NGC 2110 

55t« 

17+7 
1 ' -6 

9ll 

o 7+0.2 

1 -0.3 

14611 

4011 

0-88±S;;S 

43.311° 

0 O c >8+ 0 ' 003 
U.UOO — 0 002 

1 o + 0 - 5 

1,u — 0.4 

o.sll.t 

MCG -5-23-16 

58l| 

20±l 

Yi 

9 -1 +0.2 

CO CO 
+ 1 

4811 

0 S5+ 0 - 03 

u.oo_ 0 05 

43.911.1 

0 114+°- 004 
0.003 

9 q+0.3 
Z.9_o 2 

2 n+ 0 - 3 
Z,U -0.3 

NGC 3081 

62 ±i 

nl® 

12 ±| 

2 6+ 0 ’ 2 
Z -°-0.4 

98+^ 2 

66lH 

0 96+ 0 02 
U.DO —0 03 

43.111:1 

0 043+ 0 ' 003 

U.U^iO —q Q 03 

0.4111 

n q+o.4 

u -' 4 -o.i 

NGC 3227 

57±{ 

20±\ 

i3±| 

0 o + 01 

u,u -o.o 

147+2 

6ll 

0 94+ 0 - 01 
u,y ^-0.02 

43.011.1 

0 041+ 0 003 

u.u^i_o 003 

o.sii.i 

0 7 + 0 - 1 
u - 1 -0.1 

Circinus 

64 1? 

16±l 

12±^ 

n 4+ 0 - 3 
U, ^-0.2 

35+2 

37t S 7 

0 96+ 0 02 
u.uo_o 03 

43.411:1 

0 062+ 0 ’ 004 
u.udZ-q 003 

1 0+°' 2 
±,U -0.2 

0 9+ 0 ' 2 
u - y -o.i 

NGC 5506 

48t 6 3 

1611 

10t 2 3 

n 9+0 • 3 
U,Z -0.2 

79^5 

32+i 

0 79 + 0 - 04 
u - « y -0.05 

44 -0to.i 

0 1 30+° ’ 008 

U. loU — 0 007 

9 i+0.5 
Z>i -0.5 

1 6+ 0 ' 4 
1 ,u -0.3 

IC 5063 

6i 11 

14+ 7 

12+! 

2 - 5 -lm 

101^ 

77 H 2 

0 96+ 0 ' 02 
u.uu_ 0 04 

44 -3lo.i 

n 1 09+0.009 
U. loz _ 0 00g 

2.6111 

2 T+ 1 - 4 
z, °-1.3 

NGC 7582 

53t? 

20+ 2 

1211 

0 1+ 0-1 
u,± -0.0 

79^9 

6ll 

0 90+°'° 2 
u.uu_ 0 02 

43.5+o'.i 

0 o70+°' 005 
u.u 1 u_ 0 004 

1.411.1 

1 9 + 0-2 

i,Z - 0.2 

NGC 7674 

39±g 3 

1511 

8t 3 

11+ 0 - 5 
-‘-•-‘-- 0.6 

133l? 5 

44111 

0 16+ 0 ’ 17 
u.oo_o 20 

44 -8+°' 3 

0 33Q+ 0 - 068 
u.oou_ 0 03g 

5.3+ 2 ' 3 

0 o + l-4 

°-°-i.o 

NHBLR 

NGC 1386 

56_9 

19-5 

8+ 2 

i q + 0.3 
i,C> -0.5 

37^ 

68_5 

0.877o'.ii 

42.5+q' 3 

0 Q23+ U - UU2 
u.uzo_ 0 002 


0 4+0' 1 
u - -0.1 

NGC 3281 

6811, 

1911 


0 4+ 0-2 
U, ^- 0.2 

3811 

19±? 

o.99tS;oi 

44 -2 +q'i 

0 151 +°- 010 
u.ioi_o 008 

2 -9^o.4 

9 7+O.5 
‘ -0.4 

Cen A 

50l“ 

17l? 

loll 

Q q + 0.3 
U,, 3 - 0.2 

89+ 11 

ou_is 

3811 

O-Slto.ie 

42 - 5 to'.i 

0 021+ 0 ' 002 
u.uzi_ Q 002 

0 4+ 0-1 

0 3+ 0 - 1 
U,, 3 - 0.1 

NGC 5135 

6311 

17+5 
1 ' -2 

1211 

0-4l°o. 4 s 

71±S 

l 7 H°o 

0 97+ 0 - 01 
u -^ 1 -0.04 

43.6±°; 8 

0 079+ 0 ' 007 

u.u 1 »_ 0 .006 

1 4 + 0-4 
1 ‘^ — 0.3 

1 2+ 0 ’ 4 
i,Z -0.3 

NGC 5643 


14±J 

13±| 

o-slS:i 

sell 1 

74+8 

‘^-12 

0 97+ 0 - 02 
U -D 1 -0.04 

43.0t°' 8 

0 040+°' 004 

U.U 41 U — o.o 03 

0 6+ 0 ’ 2 
u -°-o.i 

n c+ 0.2 
U>i:> - 0.1 

NGC 5728 

6611 

17 ±l 

14+* 

0 7+ 0 - 4 

u - ‘ -0.4 

48 -6 

soli 

0 99+ 0 - 01 

u.aa_ 0 01 

43.411:1 

0 ngQ+°- 004 
U.UDo_ 0 0Q4 

1 1+ 0 - 2 

4 ■ 4 — 0.1 

1 0+°' 2 
J--D- 0.1 

NGC 7172 

69+§ 

29ti 

14+| 

0 0+ 0 ' 1 
u,u -o.o 

20 ±\ 

soil 

0 99+ 0 01 
u.aa_ 0 01 

43.411:1 

0 064+°' 002 
u,uo ^- 0.002 

1 9 + 0 ' 1 

1 s+0-i 
-*-•8-0.1 


Notes.— Torus model parameters derived from the fits with BayesClumpy. Median values of each posterior distribution are listed with their 

corresponding ±lcr values around the median. 


TABLE 5 

Torus model parameters from the global posterior distributions 


AGN 

<T torus 

Y 

No 

q 

TV 

i 

Ct 

log Air 

Tin 

T out 

H 

Type 

[deg] 





[deg] 


[erg/s] 

[pc] 

[pc] 

[pc] 

All 56^22 

Type-1 19t| 
HBLR 56+ 8 
NHBLR 64i 2 x 

1811 

19lio 

1+1 

1811 

1211 

12+1 

nil 

1311 

n c+0.0 

U,t) _o.5 

0.711:1 

o-slli 

o.4ll:l 

sill! 

H3lll 

98lH 

43lH 

45+1 

5+48 

43l®2 

48111 

0.8811:11 

o.isll.ll 

0.88ll:“ 

0.9611+1 

43.3111 

43.911:1 

43.711'.? 

43.211:1 

0 066+ U U48 
u.uoO-q 030 

0 144+ 0 - 048 
u, - l ^-0.090 

0 O66+ 0 054 
u.uoo _ 0 024 

0 O6O+ 0 006 
u.uou _ 0 042 

1 2+ u ' 4 
4 - Z -0.8 

1 6+ 0 ' 4 
4,O -0.4 

1 9+O.8 
4 - Z -0.4 

o . sil :! 

1 o + u ' 3 
4 - U - 0.6 

1 o +°'° 
4 - U - 0.6 

1 o +°' 6 

4 - U -0.3 

1 o +°' 3 
1 - u - 0.6 


Notes.— Torus parameters from the global posterior distributions of each subgroup. 
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further discussion on the Q-band excess. 

The SEDs of NGC 5506 and NGC 7172 show NIR ex¬ 
cesses compared to the model spectra as shown in Fig¬ 


ure 1 . NIR interferometric observations (Kish imoto et al, 
20091 2011 ) and NIR reverbera tion mapping (|Kishimoto 


et al. 2007 Koshida et al.||2009|) of Seyfert galaxies sug- 
gest that the AGIN torus has much smaller sublimation 
radius than expected from Eq. (1). |Kawaguchi fc Mori 


(2011) also showed in their model that if they add the rim 


darkening effect of the accretion disk, the torus inner ra¬ 
dius naturally connects to the outer disk. Thus, the NIR 
excess is readily accounted for by this c onnection. This 
is also shown in Stalevski et al. (2012) when they ap¬ 
ply the rim darkening effect to produce the model SEDs. 
However, in this study, we did not include the SED of 
a hot dust component or apply such torus geometry in¬ 
cluding rim darkening effect for the fitting as this hot 
dust remains rather unconstrained. Further discussion 
on the possi ble origins of NIR excesses can be found in 
Alonso-Herrero et al. (2011). 


4.2. General Torus Properties for Whole Sample 

In this section, we describe how we obtained the global 
distribution of the torus model parameters for the whole 
sample and the three subgroups considered here. To take 
full advantage of the data employed here, we apply a hier¬ 
archical Bayesian approach to derive information about 
the global distribution of the CLUMPY parameters for 
a given subgroup. We use a generalized beta distribution 
as the prior for each parameter (given that they are de¬ 
fined in closed intervals) and learn the hyperparameters 


of the prior using importance sampling (e.g. Brewer & 


' p _ . 

Elliott] 2014). This allows us to derive the posterior dis 


tribution for each parameter taking into account all the 
observed data that belongs to an AGN subgroup. 

In Table 5 we report the median parameters of the 
global posterior distributions for the whole sample, as 
well as for each subgroup. Note that the global posterior 
distribution of each parameter can be derived only in the 
circumstance that the prior distribution is the same for 
all the sources. This is not the case for the inclination 
angle because we use the constrained prior distributions 
for some sources as described in Section 3. For the in¬ 
clination angle, we derive the median parameters of the 
individual galaxy fits within the each subgroup. 

B ased o n interferometric observations, Kishimo to et| 
al. (2011) reported a typical torus half-light radius of 
1 pc for local AGN at MIR wavelengths. Our derived 
torus outer radius r out = 1.2+Q g pc is consistent with 
the interferometry results, although a little larger. This 
is understood because here we are considering colder dust 
within the torus than that traced by the MIR interfer¬ 
ometry, as we include data at A > 20 /am in our fits. 
These colder clumps will be generally located at larger 
radii, which explain the value of r ou t obtained from the 
global posterior distribution. The torus outer radius in¬ 
cluding cooler dust would be larger still than the value 
of Amt ^ 1.2 pc. Further studies including Herschel , 
SOFIA, and/or ALMA observations will help constrain¬ 
ing the extent of the FIR-emitting dust. Indeed, some 
pilot studies alread y showed th e importance of F IR data 


to trace cooler dust (|Ramos Almeida et al.|2011b||Garcfa- 
Burillo et al. 2014). 


4.3. Distribution of Torus Model Parameters 

The key science direction of this paper is to investigate 
the torus morphology quantitatively by deriving torus 
parameters for each subgroup and comparing them. The 
median values of the model parameters fitted to our nu¬ 
clear IR SEDs are reported in Table 4. The median 
values of the global posterior distributions of the free 
parameters cr, Y, No , g, and ry for each subgroup are 
compiled in Table 5. As discussed in Section 4.2, for the 
inclination angle i we derive the median values from the 
individual galaxy fits for each subgroup. 

Figure 2 shows the global posterior distributions of 
each physical parameter for the different subgroups. 
Black, blue, red, and green histograms show the pa¬ 
rameter distribution of all, type-1, HBLR, and NHBLR 
sources, respectively. 

As shown in Figure 2, for some of the parameters, the 
global distributions are similar for the three subgroups, 
but others are clearly different. In order to quantify t hese 
differences, we follow the same approach as in |Ramos| 
Almeida et al. (2011a ). They used the Ku llback-Leibler 
divergence (KLD^ Fullback fc Leibler|1951 ) to show that 


the joint posterior distributions of type-1 and type-2 
AGN were quantitatively different. The KLD takes into 
account the full shape of two posterior distributions to 
compare them. When the two distributions are iden¬ 
tical, the value is KLD = 0 and, the larger t he K LD 
value, the more differ ent the two distributions. |Ramos 


Almeida et al. (2011a) concluded that if KLD > 1.0, the 


two posteriors can be considered significantly different 
We calculate the KLD values for the global distributions 
of each torus parameter among the three groups. The 
values are reported in Table 6. 

We find significant differences for the parameters cr, 
y, and No between type-1 AGN and HBLR AGN. The 
differences in these parameters be tween t ype-1 a nd type- 
2 AGN were already reported in Ramos Almeida et al. 


(2011a) with larger significance, but based on fits to 


NIR and MIR photometry only, where spectroscopic data 
were not included. Besides, they did not consider in¬ 
formation from spectropolarimetry data, as we are do¬ 
ing here. Ther e fore, w e confirm the results of [Ramos 
Almeida et al.| (|2011a|) after including N band spec¬ 
troscopy to the IR photometry, which is crucial to con- 


strain the six t oru s parameters ([Alonso-Herrero et al. 
2011; Ramos Almeida et al.|2014) 
paxameterT^TT 1 


We also find that the 
and No of HBLR and NHBLR AGN 
are significantly different. Considering the average val¬ 
ues in Table 5, the tori of NHBLR AGN have larger cr, 
larger F, and larger No than those of HBLR AGN. 

There are various possible interpretations for the dif¬ 
ference in cr among the subgroups. One possibility is 
that the smaller cr could be du e to l arg er AGN lumi¬ 


nosities (receding torus model; Lawrence 


aL|2013 ). However, as shown in Table 5, the differences 


[Ri 

.tier 


1991 Ricci et 


ofmedian bolometric AGN luminosities among the sub¬ 
groups are within the uncertainties. Therefore, we con¬ 
sider that the effect of the AGN luminosity is negligible in 
this study. Another possible interpretation is a s electio n 


bias in the opt i cal typ e-1 /t ype-2 AGN se lection. Ramos 
Almeida et al.| (|20lia|) and|Elitzur| (|2012|) discussed that 
AGN classification would depend on the distribution of 


epend on the distribution of 
the obscuring material; type-1 AGN would be preferen- 
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TABLE 6 

Results of KLD test for each parameter among each subgroup 


AGN Type 

<T torus 

Y 

N 0 q 

TV 

Type-1 vs HBLR 

3.13 

1.26 

0.10 0.42 

0.10 

Type-1 vs NHBLR 

2.22 

0.22 

0.05 0.20 

4.83 

HBLR vs NHBLR 

1.78 

1.63 

0.24 0.53 

1.56 


Notes.— KLD is calculated for the global posterior distribution of each parameter among the subgroups. Values larger than 1 are shown in bold. 



Fig. 2.— Histograms of each physical parameter discussed in Section 4.3. Top panel of each figure represents the histogram of the whole 
sample. Blue/red/green filled color represents the histogram of Type-l/HBLR/NHBLR, respectively. 



Fig. 3. — Same as in Figure 1, but for the r ou t, H , and Ct parameters. 


tially selected from lower-obscuration AGN, while type- 
2 AGN (HBLR and NHBLR) from higher-obscuration 
AGN. This could be partly producing the differences in 
a that we found for type-1 and type-2 AGN, but correct¬ 
ing this effect quantitatively is extremely difficult and 
beyond the scope of this paper. 

4.4. Distribution of Covering Factor 

As shown in Section 3.1., we can derive physical pa¬ 
rameters of the torus model by combining the model pa¬ 


rameters. The individual values of these physical param¬ 
eters are reported in Table 4, and those obtained from 
the global posterior distribution for each subgroup are 
shown in Figure 3 and Table 5. 

An interesting comparison can be made between the 
geometrical covering factor of the torus model (Ct; de¬ 
scribed in equation 6) of the different subgroups and 
the average column densities derived from X-ray data 
(TVh)- NHBLRs have the largest column densities, with 
an average value of log A^h ~ 24.0 cm -2 (i.e. Compton- 
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thick), followed by HBLRs, with logTVn ~ 23.4 cm -2 , 
and type-ls, with log TVh ^ 21.8 cm -2 . Based on hard X- 
ray (50-200 keV) o bservatio ns of nearby AGN obtained 
with INTEGRAL , |Ricci et al. (2011) reported differ¬ 
ences between the X-ray reflection component of type-1 
and type-2 AGN. Type-1 and “lightly obscured” AGN 
with TVh < 10 23 cm -2 have the same X-ray reflection 
component, with reflection amplitude R ~ 0.4. On the 
other hand, “mildly obscured” AGN (TVh > 10 23 cm -2 ) 
show a clearly stronger X-ray reflection component with 
R ~ 2.2, suggesting that the central engine of “mildly 
obscured” AGN would be covered by an X-ray ref lec¬ 
ti on w all. Our results are in good agreement with |Ricci| 


et al. (2011) if we consider the Ct and TVh values for 
each subgroup. The type-1 AGN in our sample fall in 
the “lightly obscured” AGN in their study, and indeed 
they show small covering factors (Ct ~ 0.18), suggesting 
small torus X-ray reflection solid angle. The HBLR and 
NHBLR AGN subgroups would fall in the “mildly ob¬ 
scured” AGN category, and we found large covering fac¬ 
tors for them (Ct ^ 0.88 and 0.96 respectively), suggest¬ 
ing a lar ger X-ray reflection component (see also [Ricci et 


[ah][2014]). Figure 4 shows a schematic illustration of the 


torus geometrical differences among type-1 (top), HBLR 
(middle), and NHBLR (bottom). 

4.5. Torus Model Morphological Differences between 
HBLR and NHBLR 

In this section we focus on the differences between the 
modeled tori of HBLR and NHBLR. In the case of HBLR, 
we obtain smaller a values than for NHBLR, which is 
equivalent to larger torus opening angles (90° — a ). This 
implies that HBLR objects can have a larger scattering 
region, (see middle panel of Figure 4). The scattering 
region (shown schematically as a filled green bar) can 
be larger due to the larger opening angle of the torus, 
allowing more photons to be scattered, and hence polar¬ 
ized, from the BLR. We note that the larger opening of 
the ionization cone, if unresolved, will produce a slightly 
lower degree of polarization due to (partial) cancellation 
of those polarization vectors at the edges of the scatter¬ 
ing region. However, in our case, the increased amount 
of scattered photons will significantly increase the polar¬ 
ized flux, with only a very small reduction in the degree 
of polarization. To confirm this effect, we produced a toy 
polarization model assuming that the scattering region is 
a two dimensional biconical structure centered to the cen¬ 
tral engine. Then, we measure the degree of polarization 
and polarized flux, and found that the measured polar¬ 
ized flux is larger for HBLR than those for NHB LR, sup¬ 
porting our results. Miller & Goodrich (1990) also find 
the same results using a 3 dimensional cone and more 
sophisticated modeling, and more complex assumptions. 

In the case of the NHBLR, we obtain larger a val¬ 
ues than for HBLR. This means that the probability 
that scattered radiation from the BLR can be blocked is 
higher than for HBLR objects. This is also in agreement 
with the larger value of log TVh ^ 24.0 cm -2 estimated 
from X-ray observations of the NHBLR objects in our 
sample. 

To summarize, as shown in the bottom panel of Fig¬ 
ure 4, the chance to observe scattered (polarized) flux 
from the BLR is reduced by the double effect of (a) less 
scattering of the flux from the BLR (due to the reduced 





Fig. 4. — Schematic illustration of the torus geometry for type-1 
AGN (top), HBLR (middle), and NHBLR (bottom). The differ¬ 
ence in color intensity between the two bottom panels shows the 
difference in optical depth of the clumps ry, where darker color is 
larger ry. The orange region represents the BLR. The green area 
represents the media where some of the incoming BLR emission 
is scattered and then polarized (in all diagrams we describe only 
scattering from the polar scattering region, and ignore the inner 
equatorial scattering region, dominantly responsible for the polar¬ 
ized flux in Type 1 objects). The blue solid arrows represent the 
path of BLR photons and the blue dashed arrows represent the 
path of the polarized BLR photons. The observer is assumed to be 
on the left side of the torus with an inclination angle of 52°, 43°, 
and 48°, respectively (see Table 5). The only photons scattered 
along these lines of sight are shown. 


scattering area) and (b) more obscuration between the 
observer and the scattering region. Therefore, the classi¬ 
fication of an AGN as HBLR or NHBLR is probabilistic, 
and it would depend on the intrinsic properties of the 
torus, in particular of a. This could be a reasonable ex- 
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Fig. 5. — Plot comparing the / 60//25 colors measured from IRAS 
photometry and extrapolated from the clumpy torus model fits. 
Blue circle/red square/green diamond represents type-1, HBLR, 
and NHBLR AGN, respectively. The larger symbol represents the 
average value of each subgroup. The horizontal line at / 60//25 = 4 
is t he boundary between HBLR and NHBLR AGN according to 
iHeisler et al.l (fl99T ), where HBLR AGN have feo/f 25 < 4 and 
IN HBLR AGN jf 25 >4. 


planation for the lack of a hidden (polarized) BLR in 
some type-2 object^] 

However, we note that the classification of the galaxies 
as HBLR and NHBLR is mainly based on spectropolari- 
metric observations from 3-4 m telescopes, and some of 
the NHBLR could be then misclassified. Therefore, fur¬ 
ther higher sensitivity spectropolarimetry observations 
of NHBLR AGN with 8 m class telescopes such as Sub- 
aru/FOCAS and/or VLT/FORS2 are highly encouraged 
to search for the HBLR in those AGN (Ramos Almeida 
et al. in prep.). 


4.6. Inclination Anqle Effect on Detectability of HBLR 
in Type-2 AGN 

The dependence of HBLR detection in type-2 AGN on 
torus inclination angle is still a matter of debate. This 
idea arises from the observational trend that the IRAS 
60 fim to 25 /im color ratio (/ 60 // 25 ) of HBLR AGN is 
/ 60//25 < 4 on average while that of NHBLR AGN i s 
/ 60//25 > 4 (Heisler et al.| 1997| |Lumsden et al. 2001). 
Several authors have suggested that this trend is due to 
the inclination angle of the torus: the cooler and outer 
dust within the torus blocks the warm and inner hot dust 
for edge-on views, producing the high feo/f 25 , while the 
warm inner dust can be seen from the face-on views, 
reducing the value of feo/f25- Following this idea, type- 
2 AGN with low feo/f 25 would tend to be detected as 
HBLR AGN due to the more face-on view of the torus 
and vice versa. 

Here we can take advantage of the torus models fit¬ 
ted to our SEDs, which are available for each AGN and 


4 A similar explanation for the lack of HBLR detection in ~40% 
of type-2 objects, based on the distribution of dust within the torus 
and its inclination being not as simple as predicted by the unified 
model, was shown in the talk by C. Ramos Almeida at the Polar¬ 
ization & active galactic nuclei workshop held 2012 October 16-17 
at the Royal Observatory of Belgium. 


shown in Figure 1. The wavelength range covered by 
the models allows us to calculate feo /f^e color ratios for 
each source. We have also compiled feo/f 25 color ra¬ 
tios from IRAS for comparison, which probe much larger 
scales than our torus SEDs. Out of 21 sources, we ob¬ 
tained 17 IRAS feo/f 25 colors with good quality of fluxes 
(FQUAL = 3, which is the highest quality^ 

Figure 5 shows the relationship between the feo/f 25 
flux ratios obtained from IRAS and those from the torus 
model SEDs. The averages IRAS feo/f 25 flux ratio for 
type-1, HBLR, and NHBLR AGN are 4.1 ±2.5, 3.7±2.0 
and 6.4 ± 2.4, respectively, showing the previously men¬ 
tioned correlation for HBLR and NHBLR AGN, although 
with large error bars. However, when we compare the 
values obtained from the torus model SEDs, which ex¬ 
clude contamination from the host galaxy, we find that 
they are very similar for the three groups and smaller 
than the IRAS colors {feo/f 25 = 0.63 ±0.07, 0.80 ±0.48, 
and 0.75 ± 0.18 for type-1, HBLR, and NHBLR AGN 
respectively). 

These results show that the differences of IRAS feo // 2 s 
among the three subgroups are not produce d from the 
torus dust. A similar result was reported by [Al exan der| 
(2001), but using X-ray observations. Although the stan- 
dard deviations of the average values of the IRAS colors 
are large for the three subgroups, one possible explana¬ 
tion for the cooler IRAS colors of NHBLR in comparison 
with those of HBLR AGN could be dust emission from 
stronger starbursts in their host galaxies. This is in good 
agreement with previous results showing that highly ob¬ 
scured AGN tend to h ave higher star formation activity 
their host galaxies (Goulding et ^4.|[2Q12| Ichikawa et 


aT1 |2Q12a|b| [2014 |(Jastro et al.|j2U14|) . Therefore, large 
obscuration from the torus in IN HBLR AGN (as shown i 


(as snown m 

Figure 4) and higher star formation activity in the host 
galaxy could be somehow coupled. For examp le, using 
t hree-d imensional hydro dynamical simulations, |Wada fc 


Norman (2002) showed that starbursts and supernovae 
within the central 100 pc of host galaxies help lifting up 
the torus, suggesting that high star formation activity 
could influence the scale height of the torus. Comparing 
the nuclear and overall star formation activity of AGN 
with the torus obscuration is crucial to find out if they 


are coupled (e.j 

g., Imanishi et al. 2011 

Esquej et al.12014 

Ichikawa et al. 

20M). 


5. CONCLUSIONS 


We constructed 21 infrared torus-dominated SEDs in¬ 
cluding high spatial resolution NIR and MIR photom¬ 
etry, MIR spectroscopy, and Spitzer and Herschel FIR 
fluxes. By performing SED fitting using clumpy torus 
models and a Bayesian approach we derived torus pa¬ 
rameters such as the torus covering factor (Or), the torus 
inner and outer radius (r[ n and r out ), and the torus scale 
height (H). We divided the sample into subgroups based 
on whether or not they are optically type-1, type-2 with 
observational hidden broad line regions sign (HBLR), 
and type-2 without any observational broad line region 
signs (NHBLR). Our results are summarized as follows: 

1. Under the assumption of a clumpy distribution of 

5 See||Beichman et al. (1988} for the definition of FQUAL in the 
IRAS catalogs. False detections may be included when FQUAL < 
3. 
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the dust, we obtained a quantitative description of 
the torus geometry and intrinsic properties. We 
found that the median torus outer radius for the 
whole sample r out = 1.2 pc is consistent with the 
results from MIR interferometry observations. 

2. We found that the tori of Type-1 AGN have smaller 
cr, Y, TVh, and Ct than those of HBLR and 
NHBLR. Moreover, the tori of NHBLR are thicker 
and therefore have higher Ct than those of HBLR. 
These differences in the torus properties of HBLR 
and NHBLR AGN would make it more difficult to 
detect hidden BLR in NHBLR. 

3. Combining / 60//25 colors obtained from IRAS 
photometry and from torus model SEDs, we 
showed that the low / 60//25 measured fof HBLR 
using IRAS data are not due to a more face-on 
inclination of torus, but rather to star formation 
activity in their host galaxies. 
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